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Abstract 

Geo-neutrinos, electron anti-neutrinos produced in /3 decays of naturally occurring radioactive isotopes 
in the Earth, are a unique direct probe of our planet's interior. We report the first observation at more than 
3(7 C.L. of gco-ncutrinos, performed with the Borexino detector at Laboratori Nazionali del Gran Sasso. 
Anti-neutrinos are detected through the neutron inverse (3 decay reaction. With a 252.6 ton-yr fiducial 
exposure after all selection cuts, we detected ^■^^'^^{t.^i') geo-neutrino events, with errors corresponding 
to a 68.3% (99.73%) C.L. From the ln£ profile, the statistical significance of the Borexino geo-neutrino 
observation corresponds to a 99.997% C.L. 

Our measurement of the geo-neutrinos rate is 3.9^}'3(^3 2) events/(100 ton-yr). 

The observed prompt positron spectrum above 2.6 MeV is compatible with that expected from curopcan 
nuclear reactors (mean base line of approximately 1000 km). Our measurement of reactor anti-neutrinos 
excludes the non-oscillation hypothesis at 99.60% C.L. This measurement rejects the hypothesis of an active 
geo-reactor in the Earth's core with a power above 3TW at 95% C.L. 
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Geo-neutrinos (geo-Pe) are electron anti- 
neutrinos (Pg) produced in /3 decays of '^''K and 
of several nuclides in the chains of long-lived 
radioactive isotopes ^ssy ^^^^^ 232rp]^^ which are 
naturally present in the Earth. Information about 
the Earth's interior composition has insofar come 
exclusively from indirect probes: seismology only 
constrains the density profile, while geochemistry 
offers previsions based on chemical compositions 
of rocks from the upper Earth layers, chondritic 
meteorites, and the photosphere of the Sun. 
Geo-Pe's are direct messengers of the abundances 
and distribution of radioactive elements within our 
planet. By measuring their flux and spectrum it 
is possible to reveal the distribution of long-lived 
radioactivity in the Earth and to assess the radio- 
genic contribution to the total heat balance of the 
Earth. These pieces of information, in turn, are 
critical in understanding complex processes such 
as the generation of the Earth's magnetic field, 
mantle convection, and plate tectonics. 

Geo-Pe's were introduced by Eder [T] and 
Marx "Ti in the 1960's. The subject was thor- 
oughly reviewed by Krauss et al. in 1984 [3]. 
The possibility of detecting geo-neutrinos with 
large scintillator-based solar neutrino detectors was 
pointed out in Refs. [H [51 [7] . A first experimental 
indication for geo-De (~2.5 a C.L.) was reported by 
the KamLAND collaboration 

This letter reports the first observation of gco-Pg, 
performed with the Borexino detector at Laboratori 
Nazionali del Gran Sasso (LNGS). In the context, 
we also measured 9e from distant nuclear reactors 
with a mean-base line of approximately 1000 km, 
as first discussed in [W . 

Borexino is an unsegmented liquid scintillator de- 
tector built for the observation and measurement 
of low-energy solar neutrinos. The Borexino Col- 
laboration already reported the observation of ^Be 
solar-j/g [Ill[12] and the measurement of ^B solar- 

The liquid scintillator consists of 278 tons of pseu- 
documene (PC) doped with 1.5 g/1 of diphenyloxa- 
zole (PPO), confined within a thin spherical nylon 
vessel with a radius of 4.25 m. It is shielded from 
external radiation by 890 tons of liquid buffer, a 
solution of PC and 5.0 g/1 of the light quencher 
dimethylphthalate (DMP). A second spherical ny- 



^At present at James Martin Center for Nonproliferation 
Studies Monterey Institute of International Studies. 
^Also at Fermi National Accelerator Laboratory. 



Ion vessel with a 5.75 m radius segments the liq- 
uid buffer in two contiguous volumes and prevents 
diffusion of the radon emanating from the periph- 
ery of the detector close to the liquid scintillator. 
The liquid buffer is contained in a 13.7 m diameter 
stainless steel sphere (SSS). The SSS is housed in 
a 9m-radius, 16.9 m high domed water tank (WT), 
filled with ultra-high purity water, which serves as 
a passive shield against neutrons and gamma-rays. 
Scintillation light is detected by 2212 8" PMTs (the 
Inner Detector, ID). Cerenkov light radiated by 
muons passing through the water is measured by 
208 8" external PMTs (the Outer Detector, OD). 
A detailed description of the Borexino detector can 
be found in Ref. [HITi]. 

The unprecedentedly low intrinsic radioactivity 
achieved in Borexino, the high photon yield, and 
the large number of free target protons (~1.7x 10'^^) 
offer a unique opportunity for a sensitive search for 
j/e's in the MeV energy range. 

Borexino detects De via the inverse neutron /3 de- 
cay, 

Pe + P — 5^ e"*" -I- n, (1) 

with a threshold of 1.806 MeV. Some from the 
^^®U and ^-^^Th series are above threshold, while 
those from ""^K decays are below threshold. The 
positron from the inverse neutron (3 decay promptly 
comes to rest in the liquid scintillator and an- 
nihilates emitting two 511 keV 7-rays, yielding a 
prompt event, with a visible energy of £'piompt — 
Ep^ - 0.782 MeV (the scintillation light related to 
the proton recoil is highly quenched and negligi- 
ble). The free neutron emitted is typically cap- 
tured on protons with a mean time of r~256 /iS, re- 
sulting in the emission of a 2.22 MeV de-excitation 
7-ray, which provides a coincident delayed event. 
The characteristic time and spatial coincidence of 
prompt and delayed events offers a clean and un- 
mistakable signature of Vf, detection. 

In this paper we report the analysis of data col- 
lected between December 2007 and December 2009, 
corresponding to 537.2 days of live time. The fidu- 
cial exposure after cuts is 252.6 ton-yr. 

Fig. [1] shows the expected fi'prompt spectrum. It 
includes signals from geo-Ve (up to ^-^2.6 MeV) and 
reactor Df. (up to ~8MeV). For geo-i/e's we have 
used known energy spectra of (3~ decays, the chon- 
dritic Th/U mass ratio of 3.9, and the geo-Ve fluxes 
obtained from the Bulk Silicate Earth (BSE) geo- 
chemical model [16j , which predicts a detection rate 
of 2.5^0 5 events/(100 ton-yr) for geo-neutrinos in 
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Figure 1: Expected spectrum for electron anti-neutrinos in 
Borexino. The horizontal axis shows the kinetic plus the 
annihilation 1.022 MeV energy of the prompt positron event. 
Dashed line: total geo-Pe plus reactor-Pe spectrum without 
oscillations. Solid thick lines: geo— and reactor-Pe with 
oscillations. Dotted line (red): geo— Pe with the high (low) 
energy peak due to decays in the 238^ chain (SSSu-i-Zaa-pjl^ 
chains). Solid thin line: reactor-P^. See text for details. 

Borexino. 

The determination of the expected signal from 
reactor D^s required the collection of the detailed 
information on the time profiles of power and nu- 
clear fuel composition for nearby reactors. The dif- 
ferential reactor anti-neutrino spectrum, in units of 
i7e/(MeV cm^), is: 

i—1 * 

where the index r cycles over the N reactors con- 
sidered, the index m cycles over the total number 
of months M for the present data set, Tm is the live 
time in the m}^ month, is the distance of the de- 
tector from reactor r, P^m is the effective thermal 
power of reactor r in month to, the index i stands 
for the i-th spectral component in the set (^^^U, 
238u^ 239p^^ g^^^j 24ip^)^ ^j^g power fraction 

of the component i, Ei is the average anti-neutrino 
energy per fission of the component i, (j){E[,) is the 
anti-neutrino flux per fission of the component, 
and Pee is the survival probability of the reactor 
anti-neutrinos of energy Ep traveling the baseline 
Lr, for mixing parameters 9 = (Ato^jj sin022)- 
Eq. (|2]) the main contribution comes from 194 reac- 



tors in Europe, while other 245 reactors around the 
world [17] contribute only 2.5% of the total reactor 
signal. The De energy spectra, cj)i{Ej^J in Eq. ([2]), 
are taken from [TB] . Typical power fractions for the 
fuel components are: 

235^ . 238y . 239py . 241p^ _ 

(3) 

0.56 : 0.08 : 0.30 : 0.06 ^ ^ 

with a systematic error of 3.2% due to possible dif- 
ferences among the fuels of different cores and the 
unknown stage of burn-up in each reactor. For the 
thirty-five european reactors using MOX (Mixed 
OXide) technology, 30% of their thermal power was 
considered to have power fractions: 

235u . 238y . 239p^ . 241p^ _ 

(4) 

0.000 : 0.080 : 0.708 : 0.212 ^ ^ 

Information on the nominal thermal power and 
monthly load factor for each european reactor orig- 
inates from IAEA and EDF [19J. 

We use the interaction cross section Cpp for 
inverse-beta decay reaction in Eq. ([T]) from Ref. [2U] 
and the neutrino oscillations parameters {Ami2 = 
7.65 -10-5 eV^; sin^ 6*12 = 0.304) from Ref. [H] (this 
analysis includes a +0.6% contribution from matter 
effects in the approximation of constant Earth den- 
sity). The contribution of long-lived fission prod- 
ucts in the spent fuel (mainly ^''^Ru and ^''^Ce) 
amounts to 1.5% [23]. The expected reactor sig- 
nal with (without) neutrino oscillations and 100% 
detection efficiency is 5.7±0.3 events/(100 ton-yr) 
(9.9±0.5events/(100ton-yr)). The contributions to 
the estimated systematic error are summarized in 
Tab. [1] We included a 2% systematic error arising 
from the comparison between the IAEA and EDF 
reactor thermal power data. A conservative 0.4% 
systematic error on distances originates from the 
uncertainty on the Earth radius, and on the posi- 
tion of the reactor cores and the detector. 

An extensive calibration campaign with radioac- 
tive sources has been performed in Borexino. In Oc- 
tober 2008 an on-axis calibration system was used 
to place 7, /3, and a sources within the active vol- 
ume along the vertical axis. In January, February, 
June, and July 2009 other campaigns were carried 
out with an off-axis calibration system. These cam- 
paigns included AmBe, ^^Co, ^^^Ce, ^osjjg^ ssgj.^ 
54Mn, 65zn, 40k, ^OCo, and 222^.^. xhe AmBe 
source producing ^10 neutrons/s with energies up 
to 10 MeV was deployed in twenty-five different po- 
sitions allowing the study of the detector response 
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Table 1: Systematic uncertainties on the expected reactor- 
i>e signal. Sec Eq. ^ and accompanying text for details. 



Source 


Error 
[%] 


Source 


Error 

[%] 


Fuel composition 


3.2 




2.6 




2.5 


p 


2.0 


Long-lived isotopes 


1.0 




0.6 




0.4 




0.4 




0.02 






Total 






5.38 



to captured neutrons and to protons recoiling off 
neutrons. The measured light yields for gamma- 
rays following neutron capture on (2.22 MeV) 
and ^^C (4.95 MeV) at the center of the detector 
are 1060±5 photoelectrons (p.e.) ((t=42.1±0.2 p.e.) 
and 2368±20p.e. (cr=72±3p.e.), respectively. 

The stability of the detector response is contin- 
ually monitored during data taking and offline by 
means of data validation tools. Muons crossing the 
liquid scintillator can produce neutrons which are 
thermalized and captured on ^H. The mean light 
yield produced by the 2.22 MeV gamma-ray follow- 
ing neutron capture is found stable within 1%. In 
the liquid scintillator there is a measurable activity 
due to ^^°Po a decays. The mean light yield pro- 
duced by the 5.3 MeV a is measured to be stable 
within 0.5%. The stability of the overall detection 
efficiency is studied using the measured rates of cos- 
mogenic backgrounds. We clearly see the ±2% sea- 
sonal variation of the muon flux |22j . For further 
details on the detector monitoring and the methods 
of online calibrations see [T^ . 

The Geant4-based Borexino Monte Carlo (MC) 
was tuned on data from the calibration campaign. 
The expected geo and reactor D^. 's spectra shown in 
Fig. [l] were used as input to the MC code in order 
to simulate the detector response to interactions. 
The MC-generated geo-Pg and reactor-Pg spectra 
are shown in Fig. [2] where the energy is expressed 
as the total light yield (in units of p.e.) collected 
by the PMTs. 

In Borexino, the position of each event is deter- 
mined from the timing pattern of hit PMTs. The 
Fiducial Volume (FV) is determined with 3.8% un- 
certainty, based on the source calibration campaign. 
The maximal deviation from the calibration refer- 
ence positions measured at 4 m radius is 5 cm. 
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Light yield of prompt positron event [p.e.] 

Figure 2: Expected prompt positron event spectrum as ob- 
tained from the MC code, using the distribution in Fig. ^ 
as input and the selection cuts described in the text. The 
horizontal axis shows the number of p.e. detected by the 
PMTs. Primaries generated are 10^ events for both geo-Pe 
and reactor— i/e. See text for details. 

The event energy is a calibrated non-linear func- 
tion of the number of detected p.e. The total num- 
ber of p.e. collected by the PMTs, Q, depends on 
the energy, position, and nature of the events (light 
yield of 7-rays and electrons differ slightly due to 
the light quenching of low energy electrons). All 
these dependences are properly handled and well 
reproduced by the MC code, permitting us to per- 
form the analysis directly on the light yield spec- 
trum shown in Fig. [2] (rather than on the energy 
spectrum of Fig. [T]) . 

The following cuts are used for Pg's 
search: Qpi.ompt>410 p.e., where Qp 

rompt IS 

the PMTs light yield for the prompt event; 
700p.e.<(5doiaycd<l,250p.e., where Qdoiayod is 
the PMTs light yield for the delayed event; 
Ai?<lm, where Ai? is the reconstructed dis- 
tance between the prompt and the delayed event; 
20 /xs<Ai<1280 /^s, where At is the time interval 
between the prompt and the delayed event. The 
selection criterium for the reconstructed radius 
of the prompt event, i?prompt, sets a fiducial 
volume, providing a 0.25 m layer of active shielding 
against external backgrounds. The total detection 
efficiency with these cuts was determined by MC 
to be 0.85±0.01. 

An important source of background to the Pg's 
measurement is due to /3~~neutron emitters pro- 
duced by muons interacting in the scintillator, i.e. 
^Li (r=260ms) and «He (T==173ms) [24|- We re- 
ject these events by applying a 2s veto after each 
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Table 2: Estimated backgrounds for the Ue's. Upper limits 
are given at 90% C.L. 



Source 


Background 




[events/(100 ton-yr)] 


9Li-8Hc 


0.03±0.02 


Fast n's {/j's in WT) 


<0.01 


Fast n's (/z's in rock) 


<0.04 


Untagged muons 


O.OlliO.OOl 


Accidental coincidences 


0.080±0.001 


Time corr. background 


<0.026 


(7,n) 


<0.003 


Spontaneous fission in PMTs 


0.0030±0.0003 


(a,n) in scintillator 


0.014±0.001 


(a,n) in the buffer 


<0.061 


Total 


0.14±0.02 



muon crossing the liquid scintillator active volume. 
The veto inefficiency on the background from /3~- 
neutron emitters is 3 • 10^^. We tagged fifty- 
one (51) ^Li-*He candidates falling within the De 
cuts in coincidence with a positive signal from the 
above veto, corresponding to a measured ^Li-®He 
rate of 15.4events/(100ton-yr). The residual ^Li- 
*He background after the muon cut is equal to 
0.03±0.02events/(100 ton-yr). 

Fast neutrons can mimic 9^ events: recoiling pro- 
tons scattered by the neutron during its thermaliza- 
tion can fake a prompt signal, and the thermalized 
neutron capture on a proton produces a 2.22 MeV 
7-ray delayed signal. Fast neutrons contributing to 
our background can be produced by muons either 
crossing the Borexino WT or interacting in the rock 
around the detector. 

We reject more than 99.5% of fast neutrons orig- 
inated within the WT with a 2 ms veto following 
each muon crossing the WT but not the SSS. We 
have identified two (2) candidates faking a De event 
in coincidence with muons crossing the WT. Thus, 
the background from undetected muons crossing 
the WT is estimated as <0. 01 events/(100 ton-yr) 
with a 90% C.L. 

Fast neutrons originated by muons in the rocks 
surrounding the detector have an average energy 
of (i?)^90MeV and can penetrate for a few me- 
ters inside the detector, and eventually reach the 
active scintillator target. The background rate 
from these events was studied with a MC simu- 



lation, which used as input the energy spectrum 
of fast neutrons reported in Ref. 525] for the spe- 
cific case of LNGS. We estimate this background at 
<0.04events/(100 ton-yr) with a 90% C.L. 

Muons are typically identified and rejected by the 
OD, but can also be distinguished from a point -like 
scintillation event by the pulse shape analysis of the 
ID signal. Several categories of muons which might 
have gone undetected or not identified as muons, 
were studied in detail and the possibility that they 
cause a false event was considered. 

First, primary muons can mimic the prompt sig- 
nal and a muon-induced neutron a delayed sig- 
nal. This background is strongly suppressed, since 
muons depositing a visible energy below 8 MeV (our 
range of interest, see Fig. [T]), cross the WT and the 
buffer region without entering the scintillator: the 
probability that the muon-induced neutron is de- 
tected and falls within our cuts is negligible. Out 
of 0.7 millions detected buffer muons, eighteen (18) 
Pe candidates were selected, 2.5 - 10~^ per muon. 

Second, pairs of muon-induced neutrons follow- 
ing unrecognized muons can also simulate De events. 
Most muons producing more than one detected neu- 
tron have crossed the scintillator and are hence 
tagged with very high efficiency. In addition, many 
muon-induced neutrons are produced with multi- 
plicity higher than two. The requirement that, 
within a 2 ms window, every P,, candidate is nei- 
ther preceded nor followed by another event with 
neutron-like energy strongly suppresses this back- 
ground. 

The estimated combined background from 
muon-neutron and neutron-neutron coincidences is 
O.OlliO.OOl events/(100 ton-yr). 

We determined the background from acciden- 
tal coincidences by using an off-time coinci- 
dence window of 2-20 s. By simple scaling, in 
the 1260 fis wide time window used in the 
search, the number of accidental coincidences is 
0.080±0.001events/(100 ton-yr). 

To ensure the absence of any unknown time- 
correlated background we performed a detailed 
study using an off-time coincidence window of 
2 ms-2 s. We conclude that no significant time- 
correlated events are included in the data sam- 
ple, with a limit of <0.026events/(100 ton-yr) at 
90% C.L. 

Radioactivity in the nylon vessels, in the PMTs 
and, in the stainless steel can induce (7,n) and 
spontaneous fission reactions which produce MeV 
neutrons and hence mimic D/s events. On the basis 
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of the known radioactivity of the components, we 
estimate these backgrounds at <0.003 (90% C.L.) 
and 0.0030±0.0003, respectively. 

Neutrons of energies up to 7.3 MeV can also arise 
from ^^C{a,ny^O reactions following ^^°Po a de- 
cays, as investigated by KamLAND j9j. ^^"Po is 
the main contaminant in Borexino |12) . with a de- 
cay rate of 12±1 counts/(ton-day) on average for 
the present data set. Yet its abundance is many 
orders of magnitude lower than in the KamLAND 
scintillator. We estimate the probability for a ^^°Po 
a to trigger an (a,n) reaction as (5.0±0.3) x 10~®, 
from data from Ref. [26]. This source of back- 
ground, which is quantitatively very important in 
the search in KamLAND, yields an almost neg- 
ligible 0.014±0.001events/(100ton-yr) in Borexino, 
thanks to the much lower level of intrinsic ^^°Po 
background. 

Another source of background are (a,n) reactions 
due to ^^°Po decays in the buffer. ^^°Po contami- 
nation in the buffer was measured by counting a's 
from a sample of buffer fluid prepared in a vial and 
lowered at the center of the Counting Test Facil- 
ity of Borexino [27]. We obtained an upper limit 
of < 0.67mBq/kg at 90% C.L. for the contamina- 
tion from ^^'^Po in the buffer, an upper limit which 
is many orders of magnitude above the measured 
^^"Po contamination in the scintillator. A shal- 
low radial cut is very effective in removing pos- 
sible background originating from the buffer, due 
to the short distance of thermalization of the low- 
energy (a,n) neutrons. We estimate via MC simu- 
lations that the radial cut reduces 14-fold the back- 
ground induced by ^^"^Po in the buffer, correspond- 
ing to an upper limit for this background of < 0.061 
events/(100ton-yr) at 90% C.L. 

Table [2] summarizes all expected backgrounds ob- 
tained by scaling from 252.6 ton-yr fiducial expo- 
sure to 100 ton-yr for the sake of clarity. Indepen- 
dent errors are summed in quadrature. In conclu- 
sion, we expect 0.14±0.02events/(100 ton-yr) back- 
ground events in the Borexino search for P^'s. 

A total of twenty-one (21) D^s candidates pass 
all selection cuts described above. The spatial 
and time distributions of the candidates is uniform 
within the limited statistics of the observed sam- 
ple. The expected number of background events, 
in the present data set is 0.40±0.05. The signal to 
background ratio in the z/g's Borexino search is an 
unprecedented ^50:1. 

As shown by our MC (see Fig. |2]) , the light yield 
spectrum of the prompt events below 1,300 p. e. in- 
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Figure 3: Light yield spectrum for the positron prompt 
events of the 21 Pe candidates and the best-fit with Eq. ([sjl 
(sohd thick hne). The horizonal axis shows the number of 
p.e. detected by the PMTs. The small filled area on the 
lower left part of the spectrum is the background. Thin 
solid line: reactor-Pe signal from the fit. Dotted line (red): 
geo-Pe signal resulting from the fit. The darker area isolates 
the contribution of the geo-Pe in the total signal. The con- 
version from p.e. to energy is approximately 500 p.e. /MeV. 



eludes 100% of the geo-De signal and only 34.7% of 
the reactor-i^e signal. We do not know any other 
source of D^^s which could give a considerable con- 
tribution in this region. We notice that atmospheric 
and supernova relic Dg^s give a negligible contribu- 
tion. No geo-i^e are expected above 1,300 p.e. We 
note from Fig. |3] that of the total of twenty-one 
(21) candidates, fifteen (15) are in the geo-i'e en- 
ergy window below 1,300 p.e. and six (6) have a 
light yield exceeding 1,300 p.e. The 50:1 signal to 
background ratio and the clear separation of the 
two Pe sources in the light yield spectrum of the 
prompt event permit a clear identification and sep- 
aration of the number of events belonging to each 
source, and allow to establish observation of the 
geo~neutrinos, as described below. 

Experimental evidence of reactors disappear- 
ance and oscillations has been reported by the Kam- 
LAND collaboration on a mean base line of ap- 
proximately 200 km [551 [S]- Based on our cal- 
culation of anti-t/g fluxes from reactors, in the 
reactor-j^e window (Qprompt> 1,300 p.e.) we expect 
16. 3±1.1 events in absence of neutrino oscillations 
and 9.4±0.6 events in presence of neutrino oscilla- 
tions with parameters as determined in Ref. [2T]. 
The expected background in the reactor-z/g win- 
dow is 0.09±0.06. We observe in the reactor-i/e 
window six (6) events, and we conclude that our 
measurement of reactor-i/e is statistically compati- 
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Figure 4: Allowed regions for A^gco and Aroact at 68%, 90%, 
and 99.73% C.L. Vertical dashed lines: la range about the 
expected A^rcact (expected in presence of neutrino oscilla- 
tions). Horizontal dashed lines: range for Afgco predictions 
based on the BSE model in Ref. Horizontal solid red 

lines: predictions of the Maximal and Minimal Radiogenic 
Earth models. See text for details. 



ble with the expected signal in presence of neutrino 
oscillation. A statistical analysis excludes the no- 
oscillation hypothesis at 99.60% C.L. 

In the geo-De window (Qpi.ompt<l:300p.e.) we 
expect 5.0±0.3events from reactors (under the hy- 
pothesis of oscillations with the mixing parameters 
quoted above) and 0.31±0.05 background events. 
We observe in the geo-Pg window fifteen (15) can- 
didates. 

The hypothesis that the excess of events is due 
to a statistical fluctuation of the background plus 
the reactor events is rejected at the 99.95% C.L. 

Finally, we determine our best estimate of 
the geo-De and of the reactor-i/e rates with an 
unbinned maximum likelihood analysis of the 
twenty-one (21) observed candidates. Of maxi- 
mal interest is the light yield of the prompt events 
of the candidates, which, as shown above, permits 
to disentangle the two classes of events. We define 
the log-likelihood function as [30] '■ 



lnC{Nggo7 -^rcacti 5'rcacti ^Fv) — 

-^expected (-^gcoi -^rcacti '5'rcacti '5'fv)~I" 
N 

+ ^ ln[/p((5i, A^gco, iVreact, -Srcact, -Sfv) + fsiQi)] 



•^rcact 
*^rcact 



(5) 



where the index i cycles over the A^=21 candidate 
events, Qi is the light yield of the prompt event 
in p.e. for the i"^ candidate, Ngoo and Nj-eact are 
the number of geo-Pg and of reactor-i/e, the terms 
'S'rcact and S'fv account for systematic uncertainties, 
o'rcact=0.0538 and apy— 0.038 represent the frac- 
tional uncertainties on the reactor neutrino fluxes 
and on the fiducial volume described earlier in the 
text, Acxpcxtcd is the expected total number of i^e, 
fs is the spectrum of backgrounds quoted in Tab.[2j 

/p = (1 + S'fv) [/gco + (1 + Sreaci) /react] IS global 

spectrum, and /gco and /react are the individual 
spectra of the geo-v^ and of the reactor-Pg , respec- 
tively. 

Fig. [3] shows the comparison of the data with 
the best likelihood fit. Our best estimates are 
N^eo=9.9tti{1:l%^) and A^rcact=10.7t^i(li^o«) at 
68.3% C.L. (99.73% C.L.). Fig. ^ shows the al- 
lowed regions for A^geo and N^ce^f By studying the 
profile of the log-likelihood with respect to Agoo we 
have calculated that the null hypothesis for geo- 
Pe (*-e., A^goo=0) can be rejected at 99.997% C.L., 
which represents the statistical significance of the 
observation of geo-Vg reported in this paper. 

Scaling the best estimate of Agoo with 
the 252.6 ton-yr exposure, we obtain as our 
best measurement for the geo-neutrinos rate 
3.9ljj(l^;^)events/(100ton-yr). In Tab. |] we 
compare the measured rate with predictions of 
some of the most interesting geophysical models. In 
particular, we report as terms of comparison upper 
and lower bounds on the BSE models, considering 
the spread of U and Th abundances and their 
distributions allowed by this geochemical model; 
the expectation under the Minimal Radiogenic 
Earth scenario, which considers U and Th from 
only those Earth layers whose composition can be 
studied on direct rock-samples; the expectation 
under the Maximal Radiogenic Earth scenario, 
which assumes that all terrestrial heat (deduced 
from measurements of temperature gradients along 
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Table 3: Comparison the Borexino measurement of geo-i^e 
with predictions. See text for details. 



Source 




Geo-z^e Rate 
[events/(100 ton-yr)] 


Borexino 






BSE [ig 






BSE dl] 




2.5±0.2 


BSE ^ 




3.6 


Max. Radio 


jenic Earth 


3.9 


Min. Radiog 


"enic Earth 


1.6 



~20,000 drill holes spread over the World) is 
produced exclusively by radiogenic elements. 

The data presented in this letter unambiguously 
show, despite the limited statistics, the sensitivity 
of Borexino for detecting geo-neutrinos. Thanks to 
the extraordinarily low background and its unprece- 
dented 50:1 signal to background ratio obtained 
in the search, we establish observation of geo- 
neutrinos at 4.2ct C.L. The ratio between the mea- 
sured geo-Pe rate and the low-background non-Ue 
rate obtained in Borexino is ^20:1. The same ratio 
in KamLAND, as quoted from Ref. [9], is 73:276 or 
1:4, two orders of magnitude lower than in Borex- 
ino. 

The results for the geo-neutrinos rate, summa- 
rized in Tab. |3j hint at a higher rate for geo-Pg 
than current BSE predicts. However, the present 
uncertainty prevents firm conclusions. Given the 
very low background achieved in Borexino, a larger 
exposure will yield a smaller uncertainty and more 
definitive conclusions: we plan to accumulate at 
least an exposure of 1,000 tons -yr, which should re- 
sult in a reduction of the error by a factor of two. 

Finally, we investigate the hypothesis of a geo- 
reactor with a typical power of 3-10 TW at the 
Earth's core [32]. We assume an anti-neutrino spec- 
trum as detailed above, with power fractions of the 
fuel components as from Ref. [33j : 

235u : 238u ^ 0.75 : 0.25 (6) 

We set an upper bound for a 3TW geo-reactor at 
95% C.L. by comparing the number of expected 
(from reactors -I- geo-reactor and background) and 
measured events in the reactor-Pg energy window. 
Previously, this hypothesis had been studied with 
KamLAND data 0[34] obtaining a hmit of 6.2TW 
with a 90% C.L. 
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